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^ ! Abstract 

-i— > , 

The stability of a Bell-constrained half-space in compression is studied. 
To this end, the propagation of Rayleigh waves on the surface of the 

'"O | material when it is maintained in a static state of triaxial prestrain 

is considered. The prestrain is such that the free surface of the half- 
q ■ space is a principal plane of deformation. The exact secular equation is 

established for surface waves traveling in a principal direction of strain 
with attenuation along the principal direction normal to the free plane. 
As the half-space is put under increasing compressive loads, the speed 
of the wave eventually tends to zero and the bifurcation criterion, or 

Q\ [ stability equation, is reached. 

t""" Then the analysis is specialized to specific forms of strain energy 

functions and prestrain, and comparisons are made with results previ- 
ously obtained in the case of incompressible neo-Hookean or Mooney- 
Rivlin materials. It is found that these rubber-like incompressible 
materials may be compressed more than "Bell empirical model" ma- 
terials, but not as much as "Bell simple hyperelastic" materials, be- 
fore the critical stretches, solutions to the bifurcation criterion, are 
reached. In passing, some classes of incompressible materials which 
possess a relative-universal bifurcation criterion are presented. 



1 Introduction 

The works of Maurice Anthony Biot (1905-1985) cover a wide range of topics 
in mechanics and applied mathematics. Although much attention has been 
devoted to his contributions to the "acoustics, elasticity, and thermodynamics 
of porous media" pQ , his results in finite and incremental elasticity [2] were 
also far-reaching and are still relevant to many contemporary problems. For 
instance, he wrote a series of articles (summarized in his textbook [2]) on 
the surface and interfacial instability of elastic media under compression 
and his results found applications in rubber elasticity, viscoelasticity, folding 
of inhomogeneous/multilayered media, geological structures, etc. The idea 
underlying his resolution of these problems is the following: consider a media 
at rest under a finite compression; superpose an incremental inhomogeneous 
static deformation whose amplitude vanishes away from the interface; show 
that the initial compression leads to an interface deflection which is infinite; 
conclude that this condition corresponds to interface bulking or instability. 
Biot also noted that the dynamical counterparts to surface and interface 
stability analyses were Rayleigh and Stoneley wave propagation, respectively. 

This paper studies the propagation of Rayleigh waves on the surface of 
a compressed, internally constrained, hyperelastic half-space. The corre- 
sponding "bulking" or "bifurcation" criterion is derived by determining un- 
der which compressive loads the wave speed tends to zero in the secular 
equation. Biot often considered materials subject to incompressibility, an 
internal constraint which in nonlinear elasticity imposes that det V = 1 at 
all times, where V is the left stretch tensor. Here, the materials considered 
are subject to the constraint of Bell [3], tr V = 3. Both constraints are 
equivalent in infinitesimal linear elasticity (they reduce to: tr E = 0, where 
E is the infinitesimal strain tensor) but lead to quite different results at finite 
strains. In particular, the secular equation for Rayleigh surface waves cannot 
be deduced for Bell materials from the incompressible case. This equation 
was obtained in Ref. [I] as a cubic for the squared wave speed. However, this 
cubic corresponds to the rationalization of the exact secular equation and has 
spurious roots [5]; accordingly, the corresponding relevant bifurcation crite- 
rion would have to be carefully selected. Here the exact secular equation is 
found, as well as the exact bifurcation criterion. 

For incompressible materials, the Mooney-Rivlin form of the strain energy 
function Emr brings satisfactory correlation between theory and experiments 
for rubber-like materials; this function Smr is linear with respect to 7i and 
I2, the first and second invariants of the left Cauchy-Green strain tensor 
B = V 2 . For Bell constrained materials, the strain energy density for "sim- 
ple hyperelastic Bell materials" [3] is linear with respect to i<i and is, the 



second and third invariants of V. Regarding experimental results, the strain 
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energy function for "Bell empirical model" [6], Ebem = (2/3)/?o[2(3 — 22)] 4 
(where fto is a material constant) is reported as consistent with many trials 
on annealed metals such as Aluminum, Copper, or Zinc. After the equations 
governing the problem have been written and solved in Section 2 for a gen- 
eral form of the strain energy function for a Bell-constrained half-space, the 
analysis is specialized in Section 3 to the two specific forms of strain energy 
functions presented above, and the results are compared to those obtained 
by Biot for rubber-like materials. It turns out that the maximal compressive 
load that can be applied to a half-space before the bifurcation criterion is 
reached is larger (smaller) for simple hyperelastic Bell materials (Bell em- 
pirical model) than for Mooney-Rivlin incompressible materials. Also, the 
bifurcation criterion is the same for every material within each class, and an 
infinity of strain energy densities for which incompressible half-spaces admit 
such "universal" bifurcation criteria is presented in §3.41 Finally in Section 4, 
the pertinence of the notion of (in)stability for finitely deformed hyperelastic 
materials is briefly reviewed and the general interest of the Bell constraint is 
discussed, as opposed to the constraint of incompressibility. 



2 Resolution of the problem in the general 
case 

2.1 Finite pure homogeneous triaxial pre-stretch 

Let (O, xi, X2, X3) = (O, i, j, k) be a Cartesian rectangular coordinate system. 
Let the half-space X2 > be occupied by a hyperelastic Bell-constrained ma- 
terial, with strain energy density E. This material is subject to the internal 
constraint that for any deformation [31 E] , 

^ = trV = 3, (1) 

at all times, where V is the left stretch tensor. Hence, for isotropic Bell 
materials, E depends only upon i 2 and 13, the respective second and third 
invariants of V. So, E = E(z 2 ,«3), where 

z 2 = [(tr V) 2 - tr (V 2 )]/2, t 3 = det V, (2) 

and the constitutive equation giving the Cauchy stress tensor T is [3] 

T = pV + w l + w 2 V 2 , (3) 



where p is an arbitrary scalar, to be found from the equations of motion and 
the boundary conditions; and the material response functions Uq and bj<i are 
defined by 

lo = dE/di 3 , oj 2 = -is 1 dE/di 2 , (4) 

and should verify the Beatty-Hayes A-inequalities [3] 

W0OM3) < 0, W2OM3) > 0. (5) 

In the case where the material is maintained in a state of finite pure ho- 
mogeneous static deformation, with principal stretch ratios Ai,A2,A3, along 
the xi, X2,Xs, axes, the Cauchy stress tensor is the constant tensor T given 
by 

T = (j5 Ai+w + A^ 2 )i®i + (p A2 + Wo + A2W2)j®j + (p A 3 + a;o + A3W2)k(g)k;. 

(6) 
Here ujq and U2 are evaluated at 22, ^3 given by 

i 2 = A1A2 + A 2 A 3 + A3A1, i 3 = A1A2A3. (7) 

Of course, 

Ai + A 2 + A 3 = 3, (8) 

in order to satisfy ^TJ) . It is assumed that the boundary X2 = is free 
of tractions so that T 22 = 0; and that the compressive loads Pi and P3 
are applied at X\ = 00 and x% = 00 to maintain the deformation, so that 
Pi = -T on and P 3 = -T o33 . Hence, 

Po = -(wo+A!w 2 )/A 2 , P r = (A 2 -A r )(-u;o + ArA2W2)/A2, (r = 1, 3). (9) 

2.2 Incremental equations for surface waves 

Beatty and Hayes [8] wrote the general equations for small-amplitude mo- 
tions in a Bell-constrained material maintained in a static state of finite pure 
homogeneous deformation (as described in the previous subsection). These 
equations were then specialized by this author [4] to surface (Rayleigh) waves. 
The infinitesimal superposed wave is of the form 9ft{U \kx2)e lk( - Xl ~ vt >} , where 
U is an unknown decaying function. Hence the wave propagates in the direc- 
tion of the xi-axis with speed v and wave number k and is attenuated in the 
direction of the avaxis. The incremental tractions acting upon the planes 
X2 = const, are a\ x and a^ 2 , and the introduction of the scalars functions 
ti(kx 2 ) and £2(^2), defined by 

a* 21 (x u x 2 ,t) = kt 1 {kx 2 )e ik ^- vt \ a* 22 (x u x 2 ,t) = kt 2 (kx 2 )e ik ^- vt) , (10) 



allows for a compact and simple form of the equations of motion and of the 
boundary conditions. Explicitly, the equations of motion are jl] 

t[ + iX x X 2 H 2 - (AiA^C - pv 2 )U x = 0, 

t' 2 + it 1 -[b 3 (\j-\ 2 2 )-pv 2 ]U 2 = 0, 

U 2 + iX 1 X 2 1 U 1 = 0, (11) 

b 3 X 2 2 U[ + ib 3 \ 2 2 U 2 -h = 0. 

Here the prime denotes differentiation with respect to kx 2 , p is the mass 
density of the material, and 

_ —UJo + A1A2W2 

63 " A 2 (A 1 + A 2 ) >0 ' 

C al i = 2\ 2 a 5 a pu 2 - Xl(u 02 + X 2 a u 22 ) + AiA 2 A 3 (u;o3 + X 2 a u 23 ), (12) 

C = A 1 A 2 Cn + AiA 2 G 22 — Ci 2 — C 2 i — 2uq — (X 1 + X 2 )cj 2 , 

where the derivatives u r, ^>2T (r = 2,3) of the material response functions 
Uq, oj 2 are taken with respect to i? and evaluated at i 2 , i 3 given by ([7]). Note 
that the quantity 63 defined above is positive according to the A-inequalities 
(jHJ). Finally, the boundary conditions are simply 

ti(0)=t 2 (0) = 0. (13) 

2.3 Exact secular equation and exact bifurcation cri- 
terion 

The incremental Bell constraint Eq. fllip ^ suggests the introduction of a func- 
tion ip defined by 

U!(kx 2 ) = iif'ikXtX^x^, U 2 (kx 2 ) = ^(kX^x^. (14) 

With this choice, and by (TTT]) /i 1. the traction components t x and t 2 are ex- 
pressed in terms of if as: 

h = ib 3 Xl(X 1 X 2 W+^), t 2 = -b 3 XxX2<p"' J rX^X 2 {X 1 X 2 x C-b 3 X 1 X 2 -pv 2 )<p\ 

(15) 
and Eq. (rrT|) 9 reads 

b 3 X\y"" - (AiA 2 X C - 26 3 AxA 2 - pv 2 y + (b 3 XJ - pv 2 )tp = 0. (16) 

Now a law of exponential decay is chosen for if, 

cp(z) = Ae- SlZ + Be' S2 \ 3J( Si ) > 0, (17) 



for some constants A and B (it is implicit in the form of this solution that S\ 
and s 2 are distinct.) By substitution into ffTBT) . we see that the s, are roots 
of the following biquadratic, 

(bz\\)s A - (AxA^C - 2& 3 AiA 2 - pv 2 )s 2 + (M? - pv 2 ) = 0, 
s\ + s 2 2 = {XA^C - 26 3 AiA 2 - pv 2 )/(b 3 \ 2 ), s\s 2 2 = (b 3 X 2 - pv 2 )/{h\\). 

(18) 

The roots s\ and s\ of this real quadratic may be both real (and then they are 
positive because 9ft(sj) > 0) or both complex (and then they are conjugate 
because (TTBl i is a real polynomial); in both cases, sfs 2 > 0, and so by (fT8]K . 



< v < ^hXj/p. (19) 

The upper limit of this interval corresponds to the speed of a bulk shear wave 
propagating along the X\ direction. 

Now the boundary conditions (TT3]) . used in conjunction with ffTol) . and 
(EEED2, are: 



(AiA 2 " 1 s2+l)A+(A 1 A 2 T 1 s^+l) J B = 0, si(AiA 2 " 1 s^+l)A+s 2 (AiA 2 " 1 s?+l) J B = 0, 

(20) 
and the vanishing of the determinant for this linear homogeneous system of 
two equations gives the exact secular equation: 

[6s(A? - A*) - pv 2 } ^M? + (AiA^ 1 ^ - pv 2 ) ^JhXi-pv* = 0. (21) 

In the process, we used (TIBl ? ^ and dropped the factor si — s 2 . Note that 
by bringing the second term of (121 p to the right hand side and squaring, 
we obtain the cubic secular equation [1], which has spurious roots. Now for 
certain stretch ratios Ai, A 2 , A3, the speed v tends to zero in (121]) and the 
exact bifurcation criterion is deduced as 

h(\j-\ 2 2 ) + \ l \ 2 1 C = 0. (22) 

This equation defines a surface in the space of the stretch ratios which sep- 
arates a region where the homogeneous deformations of the Bell half-space 
are always stable from a region where they might be unstable. Of course, 
the critical stretch ratios must also satisfy the Bell constraint (JHJ). 

We now recast the secular equation for surface waves in a polynomial 
form for the positive quantity i], defined by [9], 

V = y/l ~ {pv 2 )/{h\\\ (23) 



as 

f(v) = ^ + V 2 + (^^ - ty - K 2 A = 0. (24) 

Clearly, at rj = (corresponding to a transverse bulk wave), we have /(0) = 
— A]~ 2 A 2 < 0; at i] = 1 (corresponding to v — 0), the secular equation tends 
to the bifurcation criterion /(l) = [6 3 (Af - A|) + AiA^ 1 C]/(6 3 A?) = 0. 

Up to this point, the setting was that of incremental surface motions 
and deformations for a general Bell-constrained half-space, maintained in 
a static state of arbitrary pure homogeneous triaxial stretch. More results 
may actually be obtained in this general setting regarding the conditions 
of existence and the uniqueness of a Rayleigh wave; this is done elsewhere 
[TO] . We now turn our attention to two specific types of Bell materials and 
compare the results obtained in plane and equibiaxial prestrains with those 
obtained for rubber-like incompressible materials. 



3 Specific forms of strain energy densities 

3.1 Simple hyperelastic Bell materials 

For simple hyperelastic Bell materials [3], the strain energy function Eshb is 
given by 

S SH B = C 1 (3-2 2 )+C 2 (l-z 3 ), (25) 

where C\ and C 2 are positive constants. The material response functions uo 
and u>2 and the quantities 63 and C provided by (JH) and (TT2]) are now 

u = -C 2 , u 2 = Ci/i 3 , h= fV?" C ' A x 3 v C = 2(C 2 + C 1 \^ 1 ). (26) 

A2{Ai + A2) 



In that context, the bifurcation criterion (1221) simplifies considerably to 

3Ai - A 2 = 0, (27) 

which is a particularly simple linear relationship between the stretch ratios 
Ai and A 2 . This bifurcation criterion is universal to the whole class of simple 
hyperelastic Bell materials because it does not depend on Ci, C 2 . This equa- 
tion delimits a plane in the stretch ratios space (Ai,A 2 ,As), which cuts the 
constraint plane (|HD along the straight segment going from the point (0,0,3) 
to the point (|, |, 0). Moreover, the analysis below shows that the region 
which is stable with respect to incremental perturbations (where there exists 
a root pv 2 > to the secular equation) is: 3Ai — A 2 > 0. In Figure 1(a), the 
plane fl27j) cuts the triangle of the possible values for the stretch ratios (jHJ) 
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(a) Region of stability (b) Surface wave speed 

Figure 1: Near-the-surface stability for 'simple hyperelastic Bell' materials. 



into two parts, of which the visible one is the region of linear surface stability 
of any simple hyperelastic Bell material. 

For the propagating surface wave, we write the secular equation ($2/^ in 
terms of 77 as: 



/ ( r/ )=^ + ^ + (l + 2Ar 1 A 2 )r / -Ar 2 A 



0. 



(2* 



As noted in the general case, /(0) < 0. At the other end of the interval 
ffT9|) . /(I) = (3Ai — A 2 )(Ai + A2) A^ 2 . So, because / is a monotone increasing 
function for rj > 0, there exists a root to the secular equation ([28]) in the 
interval [0, 1] if and only if: 3Ai — A 2 > 0; moreover, the root is unique. 

In Figure 1(b), the influence of the prestrain upon the speed of the surface 
wave is illustrated in the case of plane strain (A3 = 1). On the abscissa, Ai is 
increased from a compressive value (Ai < 1) to a tensile value (Ai > 1). The 
coordinate on the ordinate is the squared surface wave speed, scaled with 
respect to the transverse bulk wave speed, that is pv 2 / (//Af) . At Ai = 1, 
the half-space is isotropic (Ai = A 2 = A 3 = 1) because of ([H]) and the scaled 
squared speed is equal to 0.9126, the value found by Lord Rayleigh [TT] in the 
incompressible linear isotropic case. Under an increasing compressive load 
(Pi > 0, Ai < 1), the surface wave speed decreases until the critical stretch 
of (Ai) cr = 0.5 (see §3.3p . Conversely, under a tensile load (Pi < 0, Ai > 1), 
the surface wave speed increases, with the speed of the transverse bulk wave 
as an upper bound. 
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3.2 Bell's empirical model 

For Bell's empirical model materials [6], the strain energy function Ebem is 
given by 

£bem = !A)[2(3-z 2 )]1 (29) 

where /3q is a positive constant. The material response functions Uq and w 2 
and the quantities 6 3 and C provided by (TjJ and (fT2"l) are now 



Ai^2 c = r (Ai-A 2 ) 2 



Wo = , u 2 = i3 1 l3 {2(3-i 2 )}—4, b 3 = -^-, C=[2- \; o 2 i ]A x A 2 o; 2 . 



(30) 
In that context, the bifurcation criterion (1221) may be arranged as 



3 - (A ; ~ A2 | 2 - A r % = o, (si) 

4(3 -* 2 ) 1 K J 

which is the stability equation (7.10) of Beatty and Pan [12] for this problem; 
it is actually a cubic for Ai and A 2 . It delimits a curved surface in the stretch 
ratios space (Ai, A 2 , A3), which cuts the constraint plane (JSJ) into a part which 
is "unstable" (in the linearized theory) and a part which is always stable 
(in the linearized theory). This partition of the triangle of possible stretch 
ratios is visible on Figure 2(a). By extension from the plane strain case where 
A3 = 1 (treated below), we deduce that the visible part of the triangle is the 
stable one. For a clearer picture, Figure 2(b) (where the plane of the Figure 
coincides with the plane of the triangle) shows the intersection between the 
triangle and the bifurcation curve. 

For the propagating surface wave, we consider only the case where the 
underlying deformation is a plane strain such that A3 = 1. Then the Bell 
constraint (jHJ) reduces to Ai + A 2 = 2, and neither Ai nor A 2 may be greater 
than 2. Then the secular equation (l2"lj) reduces to: 



f(rj) = V * + V 2 + Ar/(2 - X 1 ) V - Af 2 (2 - A^ 2 = 0. (32) 

As noted in the general case, /(0) < 0. On the other hand, /(l) = 2(3Ai — 
2)/Af, and so, because / is a monotone increasing function for r\ > 0, there 
exists a root to the secular equation (1321) in the interval [0, 1] if and only if: 
3Ai — 2 > 0; moreover, the root is unique. 

3.3 Comparisons with incompressible rubber 

The stability of a deformed half-space made of incompressible rubber was first 
studied by Biot. He used the neo-Hookean model but noted [21 p. 165] that 
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(a) Region of stability (b) Bifurcation curve in triangle 

Figure 2: Near-the-surface stability for 'Bell empirical model' materials. 



the results were also valid for the Mooney-Rivlin model (Flavin [T3] solved 
explicitly this latter case.) Biot obtained the bifurcation criterion, showed 
that it was universal relative to both classes of materials, and computed the 
value of the critical stretch (Ai) cr at which the rubber half-space becomes 
"unstable" under compressive loads, first in the case of plane strain A 3 = 1, 
then in the case of the biaxial strain A2 = A3. Also, Green and Zerna [TH 
p. 137] found the critical stretch for the biaxial prestrain Ai = A3. In each 
case a cubic must be solved in order to evaluate the critical stretch. In the 
case of a general triaxial prestrain, the bifurcation criterion is as follows: 



Aj + A^A 2 + 3AiA 2 



A^ 



0. 



(33) 



We now show that for Bell's empirical models and for simple hyperelastic 
Bell materials, the critical stretch can be found explicitly. 

First we let the half-space made of Bell-constrained material be deformed 
in such a way that there is no extension in the 23-direction (A3 = 1). Ac- 
cording to ©2,3? this deformation is possible when the loads Pi = 2(1 — 
Ai)[Aio; 2 — ujq/{2 — Ai)] and P 3 = (1 — \i)[u>2 — ^o/(2 — Ai)] are applied at 
infinity. Then A2 = 2 — Ai by ([8]). For Bell's empirical model and for simple 
hyperelastic Bell materials, the bifurcation criterion (127|) reduces respectively 
to 

3(Ai) cr -2 = 0, and 4(Ai) cr - 2 = 0. (34) 

Then we let the half-space made of Bell-constrained material expand 
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freely in the ^-direction, so that P 3 = 0. Then we have A 2 = A 3 = (3 — Ai)/2 
by ©3 and ©, and the load P x = (3/2) (1- Ai)[AiW 2 -2u;o/(3-Ai)] must be 
applied at infinity to maintain the deformation. For Bell's empirical model 
and for simple hyperelastic Bell materials, the bifurcation criterion (|27"|) re- 
duces respectively to 

ll(Ai) cr -6 = 0, and 7(Ai) cr - 3 = 0. (35) 

Finally we consider that the Bell material is subject to a biaxial prestrain 
such that Ai = A 3 . Then A 2 = 3 — 2Ai and P x = P 3 = 3(1 — \i)[\ 1 u 2 — 
ujo/(3 — 2Ai)]. For Bell's empirical model and for simple hyperelastic Bell 
materials, the bifurcation criterion fl27|) reduces respectively to 

17(Ai) cr - 12 = 0, and 5(Ai) cr -3 = 0. (36) 

In Table 1, the numerical values for the critical stretches are given for the 
classes of Bell's empirical model (2nd column), of neo-Hookean and Mooney- 
Rivlin incompressible materials [21 [H] (3rd column), and of simple hypere- 
lastic Bell materials (4th column), in the cases of plane strain (3rd row) and 
of biaxial strain (2nd and 4th rows). It appears that rubber can be com- 
pressed more than Bell's empirical model but less than simple hyperelastic 
Bell materials, before it loses its near-the-surface stability (in the linearized 
theory) . 

Table 1: Critical stretch ratios (Ai) cr for surface instability 
Bell empirical rubber simple Bell 



Ai = A 3 0.706 0.666 0.600 
A 3 = 1 0.667 0.544 0.500 
A 2 = A 3 0.545 0.444 0.429 

3.4 A note on relative-universal bifurcation criteria 

The strain energy functions for the Mooney-Rivlin model and for the simple 
hyperelastic Bell material depend both upon two distinct material constants: 



J MR 



V 1 (I 1 - 3) + V 2 {I 2 - 3), Sshb = Ci(3 - z 2 ) + C 2 (l - z 3 ), (37) 



respectively, where T>\, T> 2 are constants and I\, 7 2 are the first two invariants 
of the left Cauchy-Green tensor B = V 2 . The fact that their bifurcation 
criteria are 'relative-universal' [15J to each class might come as a surprising 
result, but is easily understood once the strain energy functions are written 
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in terms of the principal stretches of the deformation [H] as E MR (/ 1 , J 2 ) = 
Wmr(Ai,A 2 ,A 3 ) and Eshb^,^) = Wshb(Ai, A 2 , A 3 ), where 

W MR = V x {\\ + A* + A* - 3) + V 2 (XlX 2 2 + X\\l + \l\\ - 3), 

Wshb = Ci(3 - A X A 2 + A 2 A 3 + A3A1) + C 2 (l - A!A 2 A 3 ). (38) 

Indeed, the bifurcation criterion for a general incompressible hyperelastic 
half-space may be written in terms of the first and second derivatives of its 
strain energy function W with respect to Aj (i = 1, 2) as [9], 

X 2 [W 1 + (2 - Ar 1 A 2 )W 2 ] + \\W X1 - 2A!A 2 ^ 12 + \\W 22 = 0. (39) 

When W is specialized to the Mooney-Rivlin form (I3"8"j) i. it yields the relative- 
universal bifurcation criterion (I3"3"|) . In fact, many subclasses of incompress- 
ible materials have a relative-universal bifurcation criterion. For instance, 
any incompressible material with the following strain energy function, 

W = Pi (A? + A™ + A™ - 3) + £> 2 (A?A™ + A"A™ + A"A™ - 3), (40) 

(where n — 1, 2, 3,. . . ) has the following relative-universal bifurcation crite- 
rion: 

(n - 1)A? +1 + A?A 2 + (n + l)AiA2 - A™ +1 = 0. (41) 

In particular, the bifurcation criterion ( 12"7|) . which happens to coincide with 
(I4ip when n = 1, is also valid for incompressible Varga materials (n = 1 in 

(HOD). 

Turning back to Bell-constrained materials, we note that it is a simple 
matter to write the quantities 6 3 and C in (1121) in terms of the derivatives of 
W"(Ai,A 2 ,A 3 ). We find that 

h= . ™)~ 2 W \ 2V C=(W 11 -2W 12 + W 22 )/X 3 , (42) 

^2^3^! — A 2 ) 

so that the bifurcation criterion (j22p for Bell materials is rewritten as 

W 1 -W 2 + \ 1 (W 11 -2W 12 + W 22 )=0. (43) 

When W is specialized to the simple hyperelastic Bell model ( BS]) 2 , it yields 
the relative-universal bifurcation criterion ( 12 7p . Similarly, when W is special- 
ized to the strain energy function of the Bell empirical model (T29]) . written 
as 

W BEM = f/?o[2(3 - A X A 2 + A 2 A 3 + A 3 Ai)]t, (44) 

it yields the bifurcation criterion (I3~T|) . 
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4 Concluding remarks on stability and on the 
Bell constraint 

A word of caution is needed, in conclusion, regarding the notion of instability. 
Throughout the paper, care was taken to talk of instability 'in the linearized 
theory' and of a region in the stretch ratios space where the deformed half- 
space 'might be unstable'. This is so because conclusions about the actual 
stability of a finitely deformed half-space do not necessarily come out of 
the dynamical method of surface wave analysis. As Chadwick and Jarvis 
[TB] pointed out, 'the exponential growth of a solution obtained (...) on the 
basis of a linearized theory eventually violates the assumption underlying the 
linearization'. Some authors have linked instability analysis and bifurcation 
theory, but as Guz remarked [T71 p. 268], such a comparative analysis 'is 
only qualitative and to a certain extend sketchy'. Finally, Biot [2 J adopted 
a static approach to the problem and found that at the critical stretch, the 
surface deflection (i.e. the component of the displacement normal to the 
surface) became infinite; however, it is clear that such a deformation can 
hardly be called 'incremental'. Nevertheless, it is comforting to remark that 
each approach yields the same result for the critical stretches, and to know 
that, to some extent, concording experimental results exist |18j . On the other 
hand, the linear theory of stability has its limits. As kindly pointed out by 
a referee, Chadwick and Jarvis [TB] go on to say that in a situation where 
exponential growth occurs, 'further enquiry is need to discover whether or 
not the terms initially neglected cause the motion to be stabilized.' Such a 
line of inquiry has indeed been followed since (see for instance Fu [19] or Fu 
and Rogerson [2U]), with the result that in certain cases, terms neglected in 
the linear theory do indeed contribute to a greater stability. 

Finally, the motivation for the use of the Bell constraint is now exposed. 
Using his countless experiments, James F. Bell produced a large literature 
corroborating the existence for certain metals of the constraint that now 
bears his name. A thorough background and detailed account of his research 
is given in a recent review by Beatty J2JJ. Note however that the actual 
existence of 'Bell materials' is controversial and that Bell's results have been 
criticized [221123]. From a theoretical point of view, there is a justification in 
studying classical problems (such as the one presented here) for a material 
subject to an internal constraint other than incompressibility. Indeed incom- 
pressibility is an exceptional isotropic constraint because the corresponding 
reaction stress is spherical. This property distinguishes incompressibility 
among generic isotropic constraints. For instance Pucci and Saccomandi [21] 
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proved that the deformation 

X\ = AX\ + sin XX 2 , x 2 = DX 2 , x 3 = AX 3 — cos AX 2 , 

where A, D, A are constants, is universal for all isotropically constrained 
materials, apart from incompressible materials. Also, it is always possible to 
subject a constrained material successively to a triaxial stretch followed by 
a simple shear, 

X\ = Ai-Xi + k\ 2 X 2 , x 2 = X 2 X 2 , x 3 = A3 A3, 

where the Aj and k are constants, in such a way that the following relations 
for the Cauchy stress components are universal, U\\ = cr 22 and <J\ 2 = 0, 
except when the reaction stress is spherical (these and related points are 
developed in detail by Saccomandi in Refs.[T5] or [22].) Therefore, a better 
understanding of the (theoretical) behaviour of Bell constrained materials 
gives us a better understanding of isotropic constraints and of their general 
- and not exceptional - mechanical properties. 
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List of Figures and Tables 

Figure 1: Near-the-surface stability for 'simple hyperelastic Bell' 
materials. 

Figure 1(a): Region of stability. 

Legend on graduated axes: "Ai", "A2", and "A3". 

Figure 1(b): Surface wave speed. 

Legend on graduated horizontal axis: "stretch ratio" . 
Legend on graduated vertical axis: "scaled squared speed" . 

Figure 2: Near-the-surface stability for 'Bell empirical model' ma- 
terials. 

Figure 2(a): Region of stability. 

Legend on graduated axes: "Ai", "A2", and "A3". 

Figure 2(b): Bifurcation curve in triangle. 

Table 1: Critical stretch ratios (Ai) cr for surface instability. 
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